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Investigation of spectra from the Sloan Digital Sky Survey reveals the presence 
of a region of ionized gas of > 2° diameter centered approximately at a = lO'^ST"* 
S = -00°18' (J2000) (Galactic coordinates I = 248,6 = +48). [O III] 4959,5007 
emission is particularly strong and emission from Ha and [N II] 6548,6583 is also 
detectable over a substantial area on the sky. The combination of emission line 
ratios, the close to zero heliocentric radial velocity and the morphology of the 
structure are consistent with an identification as a very nearby planetary nebula. 
The proximity of the hot, DO white dwarf PG 1034+001 further strengthens 
this interpretation. The object is: i) the largest planetary nebula on the sky, ii) 
certainly closer than any planetary nebula other than Sh 2-216, iii) the first to 
be unambiguously associated with a DO white dwarf. A parallax distance for 
PG 1034+001 would establish whether the structure is in fact the closest, and 
one of the physically largest, planetary nebula known. 

Subject headings: ISM: planetary nebulae: individual: (Hewett 1) — Stars: white 
dwarfs: individual (PG1034+001 

1. Introduction 

The availability (Abazajian et al. 2003) of a significant fraction of the spectra and images 
from the Sloan Digital Sky Survey (SDSS) (York et al. 2000) provides a unique resource for 
the investigation of a wide variety of astrophysical phenomena. The quality and scale of the 
database are such that a number of serendipitous discoveries can be expected in the coming 
years. In this Letter we report the discovery, using SDSS spectra of unrelated objects, of the 
largest known planetary nebula (PN) on the sky. 

2. Data 

In the course of a search for SDSS spectra that show the signature of two objects, at 
different redshifts, the distinctive presence of the [O III] 4959,5007 doublet, at essentially 
rest-frame wavelength, in several adjacent spectra was noted. A more targeted search em- 
ployed a simple 41-pixel (57 A) median filter to generate a "continuum" which was then 
subtracted from the original spectrum to produce a "difference" spectrum. Emission lines 
were identified in individual and composite difference-spectra using standard matched-filter 
techniques. The search revealed the presence of [O III] 4959,5007 in more than 100 spectra 
with the fiux in the [O III] 5007 line ranging from 1.6 x 10~^^ergs~^cm~^ down to the limit of 
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detectability of ~ 8 x 10 ^^ergs ^cm ^. Surface brightnesses, per square arcsecond, can be 
obtained from the fluxes measured in the spectra by dividing by the fibre area (7.1 arcsec^). 

The detections were confined to objects in a region several degrees across centered at 
approximately a = 10''37"* 5 = -00°18' or 159.3°,-0.3° (J2000). Ha, H/5, [N II] 6548,6583 
and [S II] 6718,6732 emission lines were also present in spectra in the same region. Figure 1 
shows the spatial distribution of the spectra with detectable [O III] 4959,5007 (•), Ha (o), 
and [N II] 6583 (x). The hatched area, extending to a radius of 1° from lO'^ST™ -00°18', 
indicates a region where composite spectra, derived using groups of 25 spectra without 
individual [O III] 4959,5007 detections, show unambiguous evidence of [O III] 4959,5007 
emission. Positions of objects with SDSS spectra for which no individual detections were 
obtained are also indicated (.). 

Figure 2 shows the wavelength regions containing [O III] 4959,5007 and Ha, [N II] 6548,6583 
plus [S II] 6718,6732 for a composite, continuum-subtracted, spectrum made using all the 
galaxy and quasar spectra within 0.5° of 10''37'" -00°18' (upper panels). Another compos- 
ite, using objects from a more distant arc-shaped region to the south-west is also shown 
(lower panels). The arc-shaped region is defined by radial distance 0.7° — 1.3° from 10'*37'" 
— 00°18' and angular extent 180° < PA < 315° relative to the same position. The com- 
posite from the central region displays very strong [O III] 4959,5007 and, while Ha and 
[N II] 6548,6583 are clearly visible, shows the relative weakness of the hydrogen lines. 
The composite spectrum from large radii to the south-west illustrates the fall-off in the 
strength of the [O III] 4959,5007 emission with radius and the marked variation in the 
Ha/[N II] 6548,6583 ratio with position. The spectra of stars were not included in the 
generation of the composite spectra. To avoid contamination of the composite spectrum 
from unrelated emission and absorption features, wavelength regions associated with strong 
emission and absorption features in the rest-frames of the galaxy and quasar spectra were 
excluded from the construction of the composite spectra. 

Examination of the sky-subtracted sky spectra in SDSS plates 273 and 274, which 
contain the region, show no evidence for absorption at the location of any of the emission 
hues. The lack of absorption confirms that the emission hue fluxes associated with the ~ 2° 
region, centered on 10'*37"* — 00°18', are not affected significantly by more extended emission 
on scales of ~ 5°. 

The resolution of the SDSS spectra (~ 3.lA) precludes a rehable determination of the 
radial velocity. The centroids of the [O III] 4959,5007, Ha and [N II] 6583 emission hues 

in a composite of the 38 spectra showing the strongest [O III] 4959,5007 emission give a 
heliocentric radial velocity of —5 ± 5kms^^. However, the amplitude is comparable to the 
wavelength accuracy of the SDSS spectra and variations of several tens of kms~^ are evident 



-4- 



Table 1: Emission Line Fluxes in the < 0.5° Composite Spectrum 



Species 


Wavelength 


Flux 




(iViigsiromtij 


(lU ergs ^cm 


[Ne III] 


3870 


6.5 ±1.0 


H/? 


4861 


1.0 ±0.8 


[0 III] 


4959 


10.6 ± 1.5 


[0 III] 


5007 


38.6 ± 3.0 


[N II] 


6548 


1.8 ±0.5 


Ha 


6563 


7.7 ±1.2 


[N II] 


6583 


5.1 ± 1.0 


[S II] 


6717 


1.7 ±0.7 



from spectrum to spectrum. In summary, the heliocentric velocity of the gas is consistent 
with a value of ± 20 km s~^. 

The large angular extent and small radial velocity suggest a relatively local origin for 

the ionized gas. The Galactic coordinates, / = 248, b = +48, further suggest that the gas is 
within a few hundred parsecs if the object lies within the Galactic disc. The weakness of the 
hydrogen lines and the lack of any bright early-type stars in the vicinity of the nebula rule 
out identification as an H II region. 

3. Further Investigations 

Narrow-band imaging of the [O III] 4959,5007 hues (Ac = 5008A, AA = lOOA) and 
the Ha+[N II] 6548,6583 lines (Ac = 6568A, AA = 95A) was undertaken using the Wide 
Field Gamcra (WFG) on the Isaac Newton Telescope (INT) on the nights of 2003 May 1 
and 21-27. Exposure times of 900 - 1200 s (Ha) and 3 x 900 s [O III] were used, with 
shorter exposures of 600s for companion broad-band images in g and r passbands. An area 
of roughly a square degree was imaged in both [O III] and Ha. After processing through the 
INT WFC pipeline (Irwin & Lewis 2001) there was a clear detection of complex [O III] and 
Ha nebulosity extending over the whole region, confirming the reality of the spectroscopic 
detections. The resulting continuum-subtracted, stacked [O III] and Ha images are shown in 
Figure 3. The area visible in the images is indicated by the dashed outline in Figure 1. The 
centrally concentrated distribution of [O III] 4959,5007 emission is particularly striking. A 
distinctive feature present in the [O III] image is the arc-like structure visible at center-right. 
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In the < 0.5° composite (Figure 2a,b) the [O III] A5007/A4959 ratio of 3.6 ± 0.4 and 

the [N II] A6583/A6548 ratio of 2.8 ± 0.3 are both consistent with their theoretical values, 
which depend only on atomic physics for the conditions pertaining in PN, suggesting the 
spectrum provides useful diagnostic information (although new spectra of the high surface 
brightness features are needed). H/3, which is clearly present in individual spectra with 
the strongest emission line fluxes, is barely detectable in the < 0.5° composite spectrum. 
The Ra/HP ratio of 7.5 ± 3.5 indicates some reddening (formally Ay = 0.91°!). Relatively 
weak underlying absorption in the Balmer lines would help explain the low H/3 flux while 
making little difference to the strength of Hck. The nearly reddening independent ratios of 
[N II] A6583/Ha = 0.8 ± 0.2 and [N II] A6583/[S II] A6716 = 2.9 ± 1 are both within the 
ranges observed in nitrogen-rich PN. The [O III] A5007/H/3 ratio of 38 ± 10 is very high, 
although not unprecedented in PNs (and is biased high by the selection process). Further 
support for the classification as a PN comes from the detection of [Ne III] 3869. The absence 
of detectable emission from lines of He I and He II is consistent with normal abundances 
and the S/N of the composite spectrum. In summary, the emission line properties of both 
composite and individual spectra are consistent with the properties of a PN with a relatively 
hot central star. 

Indeed, the images in the atlas of ancient PNs (Tweedy & Kwitter 1996) contain some 
strikingly similar structures. The centrally concentrated morphology evident in the distri- 
bution of [O III] 4959,5007 emission is common and large variations in emission line ratio, 
including strong [N II] and [S II] towards the outer edges of old PNs, due to interactions 
between the ejecta and the ambient interstellar medium, are often seen. The most unusual 
property of the structure reported here is the angular size, which, at > 2° diameter exceeds 
that of Sh 2-216, also long considered the closest PN, with an angular size of 1.6° (Fesen, 
Blair, & Gull 1981; Tweedy & Napiwotzki 1992). 

Since there was no plausible ionizing source in the SDSS photometric catalogue, the 
APM sky survey catalogues were used for further investigation of potential ionizing sources 
of radiation. In a 1.5° region centered on 10^37™ -00°20' one bright candidate ionizing 
source stood out at 10^37^03.875' -00°08'19.59" J2000 (Epoch 1982), with UKST plate 
magnitudes and colors of R = 13.06 and Bj — R = —0.47 (see Figure 4). Using the UKST 
plate as a reference frame the refined APM POSSl (Epoch 1952) position of this source is 
10'*37"*04.047^ -00°08'20.53" J2000, giving a high proper motion of -86±5, +31±5 masyr-^ 
The color, magnitude and proper motion suggested the object was probably a hot WD at 
around RilOO pc and a search, utilizing Simbad, reveals the source is PG 1034+001, a bright, 
V = 13.23 (Landolt 1992), DO white dwarf. A very crude kinematic age of ~ 130, 000 yr, can 
be derived from the radius of the PN, ~ 1° at a distance of ~150pc, and a typical expansion 
velocity of PN, ~ 20 km s~^. The WD is currently close to the center of the region over which 
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detectable emission is present in the SDSS spectra. The WD is also in close proximity, ^15', 
to the region showing the strongest [O III] 4959,5007 emission, however, the proper motion 
vector for the WD indicates a position for the WD to the east of the present location at the 
time the PN was formed. If the arc-like feature visible in Figure 3 represents a shock, or 
a boundary, associated with material ejected at the time the PN was formed then a lower 
limit to the kinematic age for the PN can be estimated by extrapolating the motion of the 
white dwarf back to the origin of the radius of curvature of the arc. The resulting hmit on 
the age, > 50, 000 ± 10, 000 yr is plausible but an understanding of the nature of the arc-hke 
feature is necessary to validate the argument. An age of f«100,000yr is consistent with all 
the evidence but measurement of the expansion velocity and a more complete mapping of 
the morphology of the emission will provide a much improved estimate. 



4. Discussion 

Given that all observations are consistent with a PN nature for the newly discovered 
nebulosity, we designate it Hewett 1. The proximity of the DO white dwarf PG 1034+001 
(Wesemael, Green, & Liebert 1985) makes it probable that PG 1034+001 is responsible for 
the ionizing radiation. Werner, Dreizler, & Wolff (1995) derived a spectroscopic distance for 
PG 1034+001 of ISS^Igpc and confirmation of the association would make Hewett 1 one of 
the closest PNs known (Napiwotzki 2001), the 2° diameter corresponding to a physical size of 
3.5 — 7.0 pc for the likely distance range of 100 — 200 pc. Again, Sh 2-216 provides the current 
benchmark, with a trigonometric distance of 130l8 (Harris et al. 1997) and a physical size 
of ~ 3.5 pc. A parallax determination for PG 1034+001 would establish whether Hewett 1 
is even closer. 

Notwithstanding targeted searches (Werner et al. 1997), very few PNs are known to be 
associated with non-DA (hydrogen poor) white dwarfs. Hewett 1 may be the first PN to be 
discovered associated with a DO white dwarf, although PG 0108+101 (Reynolds et al. 1987) 
and PG 0109+111 (Werner et al. 1997) have candidate nebulae, neither detection is regarded 
as secure. Since DO white dwarfs are thought to evolve from PG 1159 stars (see Dreizler 
& Werner 1996, and references therein), it will be of interest to compare the properties of 
Hewett 1 both to the PNs associated with PG 1159 stars (Napiwotzki 1999) and theoretical 
models (Gorny & Tylenda 2000). The derivation of a reliable age for the PN associated with 
PG 1034+001 will be of particular interest for constraining the time scales associated with 
the late-stages of evolution of post-asymptotic giant branch stars and the origin of PG 1159 
stars and helium-rich WDs. 
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Fig. 1. — Spatial distribution of spectra with detectable [O III] 4959,5007 (•), Ho; (o), and 
[N II] 6583 (x). The hatched area indicates a region where composite spectra also show 
unambiguous evidence of [O III] 4959,5007 emission. Positions of objects with SDSS spectra 
for which no individual detections were obtained are also indicated (•). The dashed outline 
shows the area included in the narrow-^band images of Figure 3. The location of the white 
dwarf PG 1034+001 is marked by a □. 
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Fig. 2.— Wavelength regions containing [O III] 4959,5007 and Ha,[N II] 6548,6583 plus 
[S II] 6718,6732 for a composite, continuum-subtracted, spectrum made using all the galaxy 
and quasar spectra within 0.5° of 10'*37"* — 00°18' (upper panels). Another composite, using 
objects with position angles 180° < PA < 315° lying between 0.7° -1.3° from 10''37"' -00°18' 
is also shown (lower panels). The composite from the central region displays very strong 
[O III] 4959,5007 and, while Ha and [N II] 6548,6583 are clearly visible, shows the relative 
weakness of the hydrogen lines. The composite spectrum from large radii to the south-west 
illustrates the fall-off in the strength of the [O III] 4959,5007 emission with radius and the 
marked variation in the Hq!/[N II] 6548,6583 ratio with position. 



Fig. 3. — The left hand panel shows a mosaic of 6 continuum-subtracted pointings in 
Hq;+[N II] while the right panel shows the equivalent for [O III]. The images are approxi- 
mately 0.8° on a side with North to the top and East to the left. The position of the white 
dwarf PG 1034+001 is indicated by a circle in the [O III] image. Emission with complex 
structure is evident in the central regions of the images in both passbands. A well-defined 
arc, or boundary, is visible at center-right in the [O III] image. 
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Fig. 4. — Bj — R versus R color-magnitude diagram for objects within 0.75° of 10'*37'" 
— 00°18' using the APM sky catalogue magnitudes of UK Schmidt Telescope Bj and R 
plates. All objects classified as stellar are plotted (.) with the location of the white-dwarf 
PG 1034+001 highlighted (•). 



